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ABSTRACT 


The  Korean  Seismic  Research  Station  (KSRS)  short-period  array  was 
evaluated  using  data  which  cover  eight-hour  time  periods  daily  between 
May  and  June  1973.  Most  of  the  data  is  of  poor  quality  due  to  instrumental 
problems.  The  average  noise  reduction  obtained  by  beaming  was  about  2 dB 
better  than  /n  due  to  the  almost  continuous  presence  of  propagating  coherent 
noise  at  the  station.  The  50%  detection  body  wave  magnitude  threshold  was 
found  to  be  4.35  + .45  (95%  confidence  interval)  at  the  epicentral  distance 
of  60°.  The  90%  detection  threshold  was  found  to  be  5.42  + .45  magnitude 
units. 
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INTRODUCTION 


This  report  presents  the  results  of  evaluation  of  a limited  sample  of 
digital  data  from  the  Korean  Seismic  Research  Station  (KSRS).  The  data 
consists  of  45  digital  tapes,  with  one  eight-hour  time  interval  on  each 
tape,  and  covers  the  three-month  period  of  May,  June  and  July  in  1973. 

Table  1 gives  the  dates,  time  windows  and  tape  numbers  for  all  of  the  tapes. 
The  three  most  important  observations  about  the  list  of  45  tapes  are  that 
(1)  the  eight-hour  intervals  every  other  day  provide  samples  for  only  one- 
sixth  of  the  total  three-month  interval,  (2)  the  time  of  year  is  limited  to 
late  spring  and  early  summer,  and  (3)  the  time  of  day  represented  by  different 
tapes  is  highly  variable  and  provides  samples  for  all  hours  both  day  and 
night,  although  on  different  days.  This  implies  that  diurnal  noise  variation 
may  be  observable,  but  seasonal  variations  cannot  be  observed. 

Table  2 gives  the  list  of  events  from  the  NEIS  bulletin  that  fall  within 
the  time  windows  contained  on  the  45  data  tapes.  The  events  are  distributed 
around  the  world  in  the  principal  seismic  regions.  The  asterisked  events 
were  detected.  The  nine  regions  of  best  sampling  on  this  list  are  (in  order 
of  decreasing  sample  size)  Chile,  Kuriles,  Alaska,  Fijis,  Hokkaido,  New 
Guinea,  Nevada,  Mexico,  and  New  Hebrides.  The  first,  Chile,  is  represented 
by  19  events  and  the  last,  New  Hebrides,  contains  5 events. 

The  Korean  short-period  array  is  centered  about  150  km  east  of  Seoul, 
South  Korea.  It  consists  of  19  seismometers  in  the  configuration  shown  in 
Figure  1.  They  are  located  in  a pattern  of  two  rings  around  one  central 
instrument.  The  numbers  shown  are  the  tape  channel  locations  for  each 
seismometer  site,  and  it  can  be  seen  that  number  four  is  the  only  one 
located  about  halfway  between  two  rings.  The  diameter  of  the  two  rings  are 
4.25  and  9 km  . The  greatest  spacing  of  any  two  pair  of  instruments  is  10  km 
(between  Nos.  11  and  17),  and  the  smallest  spacing  is  1.5  km  (between  Nos.  4 
and  13).  Three  principal  directional  alignments  are  emphasized  in  Figure  1. 
These  are  013°  (193°),  084°  (264°),  and  135°  (315°). 
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These  events  were  detected  at  KSRS 


The  data  quality  was  highly  variable  with  about  one-fourth  of  the  tapes 
showing  as  many  as  eight  to  twelve  channels  inoperative  or  not  operating 
normally.  Another  problem  in  data  quality  was  the  inadequate  digitizing 
(quantizing)  sensitivity  for  very  small  signals  (noise  samples  in  particular). 
This  produced  an  angular  looking  waveform  plot  and  is  probably  a source  of 
significant  errors  in  power  spectrum  computation  for  noise  samples.  These 
errors  are  probably  most  significant  at  the  higher  frequencies  (greater  than 
3Hz)  and  are  probably  tolerable  at  the  frequencies  of  importance  in  tele- 
seismic  work,  less  than  2Hz. 

In  addition  to  inadequate  quantizing  levels  many  traces  also  show 
distorted  waveforms  with  jumps  between  several  "plateaus"  and  peculiar 
looking  "one-sided"  traces  with  clipping,  dropouts  or  non-linear  compression 
of  amplitudes  in  one  direction.  These  features  are  probably  associated  with 
functional  defects  in  the  on-site  data  processing  systems.  There  are  segments 
of  recordings,  however,  where  the  signals  and  background  noise  look  com- 
pletely normal.  We  estimate,  however,  that  less  than  half  of  the  data  is 
free  of  some  kind  of  system  distortion. 
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NOISE  SAMPLES 


Six  320-second  noise  samples  were  analyzed  to  determine  the  Korean 
short-period  noise  characteristics.  All  noise  samples  preceded  strong  events 
selected  by  inspecting  seismograms  to  be  certain  that  all  19  seismometers 
were  operating  normally.  All  samples  ended  about  ten  seconds  before  the 
impulsive  start  of  a clear  first  arrival  P-phase  close  to  the  predicted 
arrival  time. 

Power  spectra  of  six  noise  samples  are  shown  in  Figure  2.  On  the  top 
part  of  the  figure  the  spectra  with  the  highest  and  lowest  absolute  noise 
level  are  shown.  The  largest  difference  between  the  two  curves  is  about  4 to 
5 dB.  The  lower  part  of  the  figure  shows  the  rest  of  the  spectra  staggered 
on  the  vertical  scale  to  avoid  visual  crowding.  All  spectra  fall  of  sharply 
towards  higher  frequencies,  and  show  secondary  peaks  at  approximately  2 and  3Hz. 

Reduction  in  noise  spectra  after  beamforming  is  shown  in  Figures  3 and  4 
for  both  the  Korean  and  NORSAR  comparable  arrays.  In  the  case  of  the  NORSAR 
site  the  spectral  data  are  "Average  Subarray  Beam"  (Barnard  and  Whitelaw, 

1972,  Figure  II-2).  In  both  cases,  spectral  noise  reduction  is  about  equal 
to  the  square-root  of  the  number  of  sensors  (10  log  N for  power  spectrum). 

An  analysis  of  the  variation  of  spectral  noise  reduction  was  done  for 

the  six  noise  samples  and  several  different  azimuths,  and  velocities  as 

listed  in  Table  3.  The  results  are  summarized  in  Figure  5.  All  fifteen 

cases  listed  in  Table  3 were  considered  at  each  frequency  and  the  median, 

minimum  and  maximum  were  calculated.  It  is  clear  from  inspection  of  Figure  5 

that  beamforming  noise  reduction  produced  a median  value  within  one  dB  of 

square  root  of  N from  1.0  to  2.0Hz.  From  about  0.5  to  1.0Hz  the  median 

was  better  than  square  root  of  N by  3 to  5 dB.  The  average  improvement  over 

1/2 

N in  the  .9  to  1.5Hz  range,  which  is  of  prime  interest  in  teleseismic 

Barnard,  T.  E.  and  Whitelaw,  R.  L.,  1972,  Preliminary  evaluation  of  the 

Norwegian  short-period  array,  Extended  Array  Evaluation  Program  Report 

No.  6,  Texas  Instruments  Incorporated,  Dallas,  Texas. 
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Figure  2.  Variation  of  noise  spectrum  with  frequency  and  date. 
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NOISE  POWER  SPECTRUM  (dB  RELATIVE  TO  1(m|ir/Hz  AT  10  Hz) 


0 1 0 2 0 3 0 4 0 


FREQUENCY  (HERTZ) 

Figure  4.  NORSAR  noise  power  spectra. 
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TABLE  III 


KSRS 

DATE 

May  0 5 
May  05 
May  05 
May  05 

May  13 
May  13 
May  13 
May  13 

May  17 
May  17 
May  17 

May  19 
May  19 

May  29 

June  30 


Array  Beams  Used  in  Noise  Analysis 


VELOCITY 

(KM/SEC) 

Infinite 

14 

14 

14 


9 

14 

14 

Infinite 


Infinite 

9 

14 


Infinite 

14 


Infinite 


Infinite 
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AZIMUTH 

(DEGREES) 


P-wave  studies,  is  about  2 dB.  The  probable  explanation  for  this  effect  also 
found  by  Prahl  et  al.  (1975),  is  the  presence  of  coherent  propagating  noise 
which  is  often  present  at  coastal  stations.  As  an  independent  check  we  also 
computed  the  noise  reduction  from  estimates  of  array  cross  power  spectral 
matrices  and  found  identical  results.  The  noise  reduction  is  approximately 
the  same  for  all  velocities  appropriate  to  teleseismic  P waves.  This  is  not 
surprising  since  the  delays  involved  are  small  compared  to  the  phase  delays 
of  slowly  propagating  0.5-1.0Hz  surface  waves.  At  frequencies  lower  than 
.3Hz  the  long-period  noise  interferes  constructively  over  the  whole  array. 


Prahl,  S.  R.,  Shen,  W.  W. , Whitelaw,  R.  L.,  1975,  Preliminary  evaluation 

of  the  Korean  seismo logical  research  station  short-period  array,  Texas 
Instruments  ALEX  (01)-TR-75-05,  Alexandria,  Virginia. 
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PREFORMED  BEAM  DETECTION 


Searches  were  made  in  predicted  time  windows  for  P-phases  using  the 
seven  beams  which  are  calculated  by  the  station  processor.  The  beam  data 
channels  are  as  indicated  in  Table  4. 

As  a result  114  detections  were  made  from  a total  NE1S  list  (Table  2) 
of  209  events.  Table  5 summarizes  all  114  detections,  giving  date, 
abbreviated  geographic  region,  seismic  region,  amplitude  (in  units  based  on 
the  same  conversion  constant,  0.488  mp/count,  used  for  seismometer  amplitude) 
and  period.  For  the  purpose  of  detectability  calculation,  the  data  were 
divided  into  two  groups.  The  first  group  includes  those  events  that  were 
within  the  preformed  beam  azimuthal  windows  and  the  second  group  includes 
those  events  that  were  in  the  two  azimuth  gaps  not  covered  by  the  preformed 
beam  (080°-100°  and  185°-205°)  and  events  at  distances  greater  than  100 
degrees.  The  magnitudes  were  reduced  to  the  epicentral  distance  of  60°  and 
sorted  into  magnitude  increment  ranges.  A maximum  likelihood  procedure 
described  by  Ringdahl  (1976)  was  used  to  fit  a detection  probability  curve 
to  the  observed  probabilities  of  detection  derived  from  Table  2.  The  results 
are  given  in  Table  6 and  Figure  6.  The  classification  of  one  event  of 
magnitude  6 was  changed  from  undetected  to  detected  since  the  failure  to 
detect  this  event  was  due  to  malfunction  in  the  recording  system  and  it 
would  have  certainly  been  detected  otherwise. 

The  body-wave  detection  thresholds  together  with  their  95%  confidence 
intervals  for  on-beam  events  at  A = 60°  were  5.42  + .45  at  the  90%  level  and 
4.35  + .45  at  the  50%  level.  The  detection  thresholds  of  "all  events"  and 
off-beam  events  are  within  two  standard  deviations  of  this  value;  therefore, 
they  are  not  significantly  different.  The  50%  detection  thresholds  agree 
well  with  that  derived  by  Texas  Instruments  group  (Prahl,  Shen,  and  Whitelaw, 
1975).  The  90°  threshold  is  considerably  higher  than  at  other  arrays 


Ringdahl,  F,  1976,  On  the  estimation  of  seismic  detection  thresholds.  Bull. 


TABLE  IV 


Preformed  Beams 


Channel  # 

Azimuth 

(Degrees) 

Velocity 

(km/sec) 

32 

030 

13.0 

■ 33 

150 

13.0 

34 

240 

13.0 

35 

270 

13.0 

36 

300 

13.0 

37 

1 / 

330 

13.0 

38 

360 

13.0 

I 


TABLE  V 


Detections  from  NEIS  Bulletin 
(114  out  of  Possible  209) 


DATE 

day  event 

(Julian  (Table 
Calendar) 

GR/SR 

AMP/PER 

# OF 
GOOD 
SEIS  . 

4/29 

119 

1 

HON/ 01 

18/2.0 

19 

2 

BAN/ 01 

60/1.5 

19 

3 

PHI/02 

24/1.3 

19 

4 

PAN/ 04 

35/1.1 

19 

7 

HDK/08 

23/2.0 

19 

5/0  5 

125 

15 

SOL/02 

430/1.2 

19 

19 

IRA/06 

130/0.9 

19 

20 

HON/08 

110/0.5 

19 

23 

HOK/09 

38/1.4 

(19) 

5/13 

133 

28 

ASC/02 

38/1.5 

19 

29 

BON/05 

100/1.0 

19 

5/17 

137 

36 

KUR/06 

95/0.8 

19 

37 

SIN/10 

240/1.2 

19 

5/19 

139 

38 

MAR/ 18 

48/0.8 

19 

39 

PER/19 

310/1.2 

19 

5/21 

141 

40 

CAR/ 12 

40/0.9 

(19) 

5/25 

145 

44 

WPA/09 

90/1.5 

19 

46 

ORE/12 

35/1.0 

19 

47 

ALA/ 13 

194/0.9 

19 

5/29 

149 

51 

MAS/12 

57/1.2 

19 

52 

SSI/13 

30/0.8 

(19) 

54 

BAF/16 

22/1.0 

(19) 

55 

PBB/17 

125/1.0 

19 

5/31 

151 

56 

SUM/12 

50/0.8 

19 

58 

RAT/18 

112/0.8 

19 

6/02 

153 

59 

IRA/14 

32/0.8 

18 

6/04 

155 

60 

PHI/10 

34/0.9 

18 

61 

NBR/13 

39/1.0 

18 

62 

ECU/13 

20/0.8 

18 
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TABLE  V (Continued) 


Detections  from  NEIS  Bulletin 
(114  out  of  Possible  209) 


DATE 

, °AY  EVENT 

Calendar! Tab 

GR/SR 

AMP/PER 

# OF 
GOOD 
SEIS. 

6/06 

157 

67 

NEV/13 

2200/0.8 

7 

69 

SOL/14 

150/1.0 

14 

70 

HIN/16 

39/0.6 

14 

71 

CEL/17 

34/1.2 

19 

6/08 

159 

74 

BON/ 12 

65/0.8 

7 

75 

NGU/14 

65/0.7 

7 

6/10 

161 

80 

CEL/16 

46/1.0 

(19) 

82 

CAR/17 

47/1.3 

(19) 

83 

GUA/18 

36/1.5 

(19) 

6/12 

163 

85 

SOL/12 

150/1.0 

3 

86 

KAM/14 

570/0.8 

7 

87 

HOK/15 

386/0.8 

6 

6/14 

165 

89 

FLO/ 11 

1450/1.2 

(19) 

90 

SEV/13 

20/0.4 

(19) 

91 

SUM/15 

123/0.9 

(19) 

6/16 

167 

95 

BKL/12 

104/1.2 

14 

97 

ORE/15 

495/1.5 

14 

98 

BAN/ 15 

130/0.6 

14 

6/18 

169 

99 

HOK/13 

60/1.0 

100 

KUR/14 

40/0.9 

8 

101 

KUR/15 

66/1.2 

11 

104 

HOK/18 

315/1.0 

14 

6/20 

171 

105 

HON/ 15 

38/0.8 

(19) 

107 

KER/18 

115/0.8 

19 

6/22 

173 

108 

FIJ/17 

95/1.3 

19 

110 

NHE/18 

26/1.1 

(19) 

111 

MAR/ 19 

1807/0.8 

19 

112 

TON/ 20 

414/1.4 

19 

6/24 

175 

114 

MAR/ 17 

400/2.0 

19 

115 

KUR/18 

58/0.8 

(19) 

116 

NGU/18 

81/1.4 

(19) 

117 

KUR/19 

133/1.2 

19 

118 

KUR/20 

238/1.5 

19 

119 

NOR/21 

54/2.6 

19 

6/26 

177 

122 

KUR/18 

426/1.0 

19 

123 

HOK/18 

84/1.0 

(19) 

124 

CRE/19 

193/1.0 

19 

TABLE  V (Continued) 

Detections  from  NEIS  Bulletin 
(114  out  of  Possible  20Q) 


DATE 

DAY 

(Julian 

Calendar) 

EVENT 

(Table 

GR/SR 

AMP/PER 

# OF 
GOOD 
SEIS . 

6/28 

179 

126 

SCR/14 

412/1.0 

19 

127 

MOL/ 15 

58/1.0 

19 

128 

CHI/15 

36/2.8 

(19) 

130 

NGU/17 

75/1.1 

19 

131 

NEV/19 

39/1.0 

(19) 

133 

BAN/ 20 

138/1.2 

19 

6/30 

181 

134 

GUA/08 

84/1.2 

(19) 

7/0  2 

183 

141 

AEG/12 

22/1.0 

(19) 

142 

TAD/ 13 

38/1.0 

19 

143 

ALA/14 

32/0.9 

19 

145 

KAM/15 

25/0.8 

19 

7/04 

185 

148 

ALA/ 13 

30/0.7 

19 

150 

ICB/17 

34/1.0 

(19) 

151 

HON/ 17 

65/0.9 

19 

7/06 

187 

152 

NIC/05 

23/0.9 

19 

154 

CBB/07 

13/0.8 

(19) 

155 

CHI/04 

58/1.5 

19 

156 

CHI/11 

11/0.8 

(19) 

157 

CHI/12 

23/2.0 

19 

7/0  8 

189 

159 

OKH/15 

76/0.6 

(19) 

160 

LEE/17 

38/1.3 

19 

161 

CHI/18 

47/1.4 

19 

7/10 

191 

163 

CHI/15 

47/2.0 

19 

165 

KUR/18 

31/1.2 

19 

7/12 

193 

169 

CHI/14 

100/2.2 

19 

170 

CHI/15 

15/1.0 

(19) 

171 

CHI/16 

67/1.5 

19 

172 

CHI/16 

106/1.4 

19 

173 

PHI/19 

55/1.0 

19 

7/14 

195 

174 

GRE/13 

60/0.7 

19 

175 

TIB/14 

326/1.8 

19 

7/16 

197 

176 

FLO/13 

55/1.6 

19 

178 

MEX/18 

31/1.1 

19 

179 

TIB/20 

277/1.2 

19 

7/20 

201 

183 

CER/10 

49/1.6 

17 

7/22 

203 

188 

SSI/15 

12/1.0 

10 

189 

KER/15 

11/1.0 

10 

190 

PHI/15 

28/0.8 

(10) 

19  3 

HON/18 

23/1.2 

10 
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DATE 

7/24 


7/26 


TABLE  V (Continued) 

Detections  from  NEIS  Bulletin 
(114  out  of  Possible  209) 


DAY 

(Julian 

Calendar) 

EVENT 

(TaDle 

2) GR/SR 

AMP/PER 

# OF 
GOOD 
SEIS 

205 

197 

KUR/15 

98/1.1 

(10) 

19  8 

SOL/15 

44/1.2 

(10) 

199 

CHI/20 

264/1.3 

10 

207 

202 

MAC/14 

45/1.5 

13 

203 

AND/14 

62/1.2 

13 

205 

AND/15 

26/1.1 

13 

207 

AND/19 

91/1.2 

13 

208 

AND/20 

49/1.5 

13 

209 

AND/ 20 

401/1.2 

13 
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TABLE  VI 


Incremental  Detection  Ratios 


Mag . Range 
(Center 
Value) 

All  Events 

Events  In 
Preformed 
Beams 

Events  Outside 
Preformed 
Beams 

3.6 

.00 

.00 

— 

3.8 

.29 

. 33 

. 00 

4.0 

.43 

.75 

.00 

4.2 

.53 

.71 

.00 

4.4 

. 31 

.29 

.33 

4.6 

.58 

.54 

.71 

4.8 

.60 

.59 

.63 

5.0 

.86 

.83 

1.00 

5.2 

.88 

.83 

1.00 

5.4 

.93 

1.00 

.87 

5.6 

.91 

1.00 

.75 

5.8 

1.00 

1.00 

1.00 

6.0 

1.00* 

1.00* 

— 

6.2 

1.00 

1.00 

— 

* Changed  from  .5  to  1.0  justification  in  text 
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Figure  6.  Observed  incremental  detection  probabilities  and  the  maximu: 
likelihood  detection  probability  curve  fitted  to  data.  Filled  circles 
show  the  50%  and  90%  detection  levels. 


Figure  6 (Continued).  Observed  incremental  detection  probabilities  and  the 
maximum  likelihood  detection  probability  curve  fitted  to  data.  Filled  circles 
show  the  50%  and  90%  detection  levels. 


Figure  6 (Continued).  Observed  incremental  detection  probabilities  and  the 
maximum  likelihood  detection  probability  curve  fitted  to  data.  Filled  circ 
show  the  50%  and  90%  detection  levels. 


r 


(Dean,  1971;  Barnard  and  Whitelaw,  1972).  One  can  attribute  this  to  the 
high  noise  level  at  K.SRS  which  is  about  three  times  the  average  noise 
amplitude  at  NORSAR.  Although  data  quality  is  often  poor,  only  a few  events 
were  not  detected  because  of  instrument  malfunctions;  these  do  not  effect 
the  threshold  determination  appreciably. 


Dean,  W.  C.,  1971,  Detection  threshold  fo  the  LASA/SAAC  system,  Seismic 
Array  Analysis  Center  Report  No.  3,  Teledyne  Geotech,  Alexandria, 
Virginia. 
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UNFILTERED  SIGNAL  BEAMFORMING 


Fifteen  of  the  best  events  were  selected  for  unfiltered  signal  and  noise 
beamforming  in  the  precalculated  direction  and  azimuth  of  the  event.  In 
Table  7 we  give  the  signal  loss  at  beaming  and  the  average  wide-band  noise 
reduction  in  dB.  The  results  show  that  the  signal  loss  is  very  small 
indicating  that  the  signal  is  similar  across  the  array.  The  overall  noise 
amplitude  reduction  is  close  to  /N  (12.8  dB) . This  is  probably  because  the 
noise  cancellation  in  wide  band  use  is  dominated  by  noise  amplitudes  at 
the  spectral  peak  where  at  the  low  frequency  side  the  noise  reduction  is 
worse  than  /n  while  at  slightly  higher  frequencies  it  is  better  (see 
Figures  3 and  5),  the  overall  result  is  therefore  close  to  /n. 
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TABLE  VII 


Beamforming  to  Selected  Events 


EVENT 

AZ. 

VEL. 

SIGNAL 

LOSS 

dB 

NOISE 

REDUCTION 

dB 

SOL/12  5/02 

143.5 

15.16 

0.45 

11.48 

IRA/125/06 

288.6 

15.62 

0.85 

12.18 

BON/133/05 

130.1 

10.32 

0.39 

10.70 

KUR/1 37/06 

059.4 

8.95 

1.35 

14.37 

SIN/137/10 

290.2 

12.90 

0.31 

12.73 

PER/139/19 

047.2 

79.76 

0.57 

13.78 

ALA/145/13 

047.8 

14.70 

0.25 

12.38 

PBB/149/17 

039.9 

72.16 

1.09 

14.55 

RAT/151/18 

052.0 

13.02 

0.46 

10.10 

TON/173/20 

126.6 

20.56 

0.19 

14.90 

MAR/175/17 

137.0 

10.32 

0.33 

8.99 

KUR/175/20 

062.2 

8.53 

1.00 

10.70 

SCR/179/14 

136.1 

16.62 

0.26 

12.73 

TIB/195/14 

279.0 

12.72 

0.02 

13.26 

TIB/197/20 

278.7 

12.75 

0.71 

9.22 
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SUMMARY 

Evaluation  of  the  Korean  Seismic  Research  Station  led  to  the  following 
conclusions . 

1. )  The  average  wide  band  rms  noise  reduction  is  about  2 dB  better  than 
/N  due  to  coherent  propagating  noise  at  frequencies  below  1.5Hz.  At 
frequencies  above  1.5Hz  the  noise  is  incoherent. 

2. )  The  50%  body  wave  detection  threshold  is  at  = 4.35  + .45  at  the 
epicentral  distance  of  60°.  The  90%  detection  threshold  is  at  = 5.42 
+ .45.  The  limits  indicated  above  are  at  the  95%  confidence  level. 

3. )  Instrument  and  system  malfunctions  reduce  the  utility  and  reliabil- 
ity of  the  data  at  this  array  in  the  May  - June  1973  time  interval 
greatly. 
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